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TECHNICAL NOTE 2138 


ANALYTICAL AND EXPERIMENTAL INVESTIGATION OF ADIABATIC 
TURBULENT FLOW IN SMOOTH TUBES 
By Robert G. Deissler 


SUMMARY 

Equations were derived for the prediction of velocity distribu- 
tions for fully developed adiabatic turbulent flow In smooth tubes; 
both the incompressible- and compressible-flow cases were treated. 

The analysis produced a single equation that represents flow in both 
the conventional buffer layer and the laminar layer. By graphical 
integration of the velocity-distribution equations developed, a 
dimensionless flow-rate parameter was obtained and plotted. Use 
of the flow-rate parameter permits the prediction of pressure 
gradients along a tube for various flow rates. 

In order to check the analysis and to determine values for 
the constants appearing In the equations, tests were conducted to 
determine velocity distributions in air flowing without heat transfer 
In a smooth tube having an inside diameter of 0.87 inch and a length 
of 87 Inches. Data were obtained for fully developed velocity dis- 
tributions and for developing velocity distributions at various 
distances from the entrance of the tube. 

The results for fully developed flow were correlated by using 
conventional dimensionless velocity and distance parameters, and 
agreed closely with those of Nikuradse and other investigators. 

The plots of the equations and of the flow-rate parameter agreed 
well with the data when appropriate values of the two experimental 
constants appearing in the equations were used. 


INTRODUCTION 

Much empirical work has been done on the flow of fluids in 
tubes, and it has long been possible to predict fluid-flow pressure 
drops with a fair degree of accuracy by using experimentally deter- 
mined friction factors (reference 1, pp. 232-272). Somewhat less 
work has been dene in an effort to understand the mechanism of 
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turbulent flow in tubes and to develop a theory that is in complete 
agreement with the measurements. The Karman similarity theory, which 
is perhaps the best-known contribution, satisfactorily predicts 
the velocity distributions for fully developed incompressible flow 
except in the vicinity of the wall (reference 2). 

In an investigation made at the NACA Lewis laboratory, a new 
equation was developed that gives the relation between velocity and 
distance from the wall for both the laminar and the so-called buffer 
regions. For completeness, the Karman-Prandtl development for the 
velocity distribution at a distance from a wall (turbulent region) 
is also included. In the present analysis, this development is ex- 
tended to compressible flow, that is, to the case in which variation 
of temperature across the tube due to high velocities is appreciable. 

Because of the scarcity of velocity-distribution data, 
especially for conditions under which compressibility might be 
appreciable, investigations were made to determine velocity dis- 
tributions in a tube at various distances from the entrance and 
for fully developed flow. The results of these investigations, 
reported herein, are used to check the analysis. 


ANALYSIS 

General Turbulence Theory 

During turbulent flow through a tube, portions of the fluid 
move about in random fashion. Inasmuch as a transverse-velocity 
gradient exists, some portions enter regions of various mean axial 
velocities. Momentum is then transferred from one portion to 
another and a shear stress, in addition to the viscous shear stress. 
Is produced. 

By analogy with the law for viscous shear stress T v = p du/dy, 
the equation for the shear stress produced by turbulence is often 
given in the following form: 


T 


t 



( 1 ) 


where pc is comparable to the viscosity for viscous shear p, 

c is the coefficient of eddy diffusivity, the value for which 
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la determined by the amount and the kind of turbulent mixing at a 
point. (The symbols used in this report are defined in the appendix.) 

The analogy between shear stresses produced by viscosity arid 
those produced by turbulence is not quite exact because the mechanism 
for momentum transfer is somewhat different for the two conditions. 

The chief difference probably is that in the case of viscous shear* 
momentum transfer takes place suddenly at the instant the molecules 
collide; whereas in the case of turbulence, portions of the fluid 
can continuously transfer their momentum as they travel transversely 
inasmuch as they continuously act on one another. This difference 
can, however, be absorbed in the value of £, which is descriptive 
of the turbulence mechanism, so that equation (l) should still be 
valid. 

The total shear stress T may be obtained by adding the vis- 
cous shear stress to the turbulent shear stress as follows: 


t = n du + pg du 
dy dy 


( 2 ) 


In order to make practical use of equation (2), € must be 

evaluated for each portion of the flow. Because the actual mech- 
anism of turbulent exchange of momentum is very complicated, the 
method of dimensional analysis is used and the constants obtained 
in the analysis are evaluated by experiment. From consideration 
of the various factors on which e might depend, the following 
functional relation is assumed: 


d 2 u 


d^u 


£ = f f u V H 

V ’ P ' dy> *y 2 ' dy 3 ' 


Although the turbulence mechanism might conceivably depend on p/p 
the influence of this factor is assumed negligible and the validity 
of this assumption is experimentally checked. 

Flow at a distance from wall . - It is shown by von Karman 
(reference 2) that for flow at a distance y from a vail the 
shear stress is given by 
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( 3 ) 


where i a a constant of proportionality experimentally deter- 

mined. Viscous shear stress is neglected in the derivation of this 
equation. Combination of this equation and equation (2), with the 
viscous shear stress neglected because it is small except in the 
region close to the wall, gives the following expression for local 
eddy diffusivity for flow at a distance from the wall: 


( 4 ) 


Equation (4) could also have been obtained from the general func- 
tional relation for € by assuming that only the first and second 
derivatives of the velocity with respect to distance are important 
in determining the value of €, and by applying dimensional anal- 
ysis. This fact indicates that the eddy diffusivity, or turbulent 
transfer of momentum, at a point away from the wall is chiefly 
dependent on the velocities in the vicinity of the point relative 
to the velocity at the point and is independent of the magnitudes 
of u and y. The fact that € at a point away from the wall is 
dependent more on the velocity distribution than on the magnitude of 
the velocity u at the point may be illustrated by noting that at the 
smooth entrance to a tube, where the velocity profile is uniform over 
most of the cross section, turbulence at a point away from the wall 
is negligible compared with that farther down the tube at a point at 
which the mean axial velocity is equal to that at the entrance, but 
at which the flow profile is fully developed. This change from zero 
eddy diffusivity at the entrance to an appreciable degree of turbu- 
lence farther down the tube has been clearly shown by experiments 
with dye in a stream (reference 1, (Frontispiece.)). In the remainder 
of the analysis, equation (3) is used in calculating flow at points * • 

distant from the wall. 
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Flow In vicinity of vail . - Although the effects of the magni- 
tudes of u and y on € can be neglected in considering flow at 
points distant from the wall, it appears that they must be accounted 
for in considering flow close to the wall. The experimental data 
available indicate that the turbulent transfer of momentum and thus 
the turbulent shear stress become very small in the region close to 
the wall, so that near the wall practically all shear stress is 
produced by viscous action and the velocity u is very nearly a 
linear function of y (reference 3). The second and possibly 
higher velocity derivatives therefore approach the constant value 
zero in the vicinity of the wall and the first derivative approaches 
a constant, which is defined when u and y are given. As pre- 
viously stated, the effect of n/p on € is assumed negligible. 

As a first approximation, the functional relation for € is there- 
fore written as 


€ = f (u,y) 


The range of values of y for which this approximation is suffi- 
cient will be experimentally determined. From dimensional analysis, 


£ = n 2 uy 


(5) 


where n 2 is a constant of proportionality to be determined 
experimentally. Substitution of equation (5) into equation (2) 
gives 


T 



+ n 2 puy 


du 

dy 


( 6 ) 


In the remainder of the analysis, equation (6) is used in calculating 
velocity distributions close to the wall. 
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Equations for Incompressible Flow 

Flow In vicinity of wall . - In order to obtain equations that 
are In dimensionless form, the following commonly used dimensionless 
quantities are employed: 


u + = E_ 


Then 


,,-f - fVPO 

" W 7 


du ^ T 0 du + 
dJr “O dy + 


d£u = T 0 Vy Pp d 2 u + 
-2 2 / 2 
dy ^0 /Po dy + 


On substitution of equations (7) to (9), equation (6) becomes 


T n = — — + n 2 P„Vt u + ^ 


'0 ‘ ^0 


f + T 0 du + 




0 u O' 


°0 Vt^Oq dy + 


vnore p and |i have been replaced by and p^, respectively, 

because the density and the viscosity are the same at all points for 
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incompressible flow. The shear stress T has been approximated by 
Tq because only the region close to the wall Is being considered 
and the variation in T for this region is slight. 

By rearrangement and cancellation of like terms, the preceding 
equation simplifies to 


- 


n 2 u + y + 


= 1 


du 


( 11 ) 


This equation Is a first-order linear differential equation having 
the solution 


r \x 


y + = Ce 


/ *f* % 2 / 2 

(null (nu ) 


+ e 


-(nun 2 

e 2 du + 


U 0 


When y + = 0, then u + = 0 and 




^ 0 


f +\2 
-( nu ) - 

e 2 du + = 0 


Therefore C = 0 and 




+\2 


y + = e 2 


e 2 du + 






0 
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This equation can be rewritten as 


y T = i 
n 



r-nu 


-(sail 2 


2 d(nu + ) 


/ + \2 
l(»u ) 


1 e 2 


VS 


( 12 ) 


vhere 


/ + v< 

*(nu ) 


-L 

A/^bt 


e 


is the normal error function of nu + . (See 


n£r^ *?? yeiation between u + and y+ for incompressible 
flow near the wall is given in equation (12). 

. , at _a distan ce from wall with shear stress uniform across 

Jjetube. - For flow at a distance from the wall, the analyTiT't'hat 
was developed by von Karrnan for the region near a wall (effect of 
viscosity neglected) and applied by Prandtl to the region at a dis- 

* ' "I 1 * 18 P reaan tsd here. This application involves the 
assumption that the shear stress is constant across the tube an 
assumption that will later be shown to be sufficiently accurate for 

:L i8 ^? liCati0n - If thla a88U3I P*ion is made and expressions (9) 
and (10) are substituted into equation (3), v ' 




( 13 ) 


This equation can be integrated to give 


, *. V. 
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At the wall, the velocity gradient 1 b very large compared with that 
at a distance from the wall eo that dy + /du + can he considered 
equal to zero when y + = 0, giving = 0 and 


= i 

K , + 
4u 


(14) 


After the variables are separated, this equation can he integrated 
again to give 


u+ = ir log e y+ + c 


(15) 


where C is a constant of integration, the value of which may he 
found, when the ranges of applicability of the equations for flow 
close to a wall and at a distance from a wall are known. The ranges 
of applicability fcr the equations are to be experimentally deter- 
mined. The relation between u + and y + for incompressible flow 
at a distance from the wall when constant shear stress across the 
tube is assumed is given in equation (15). 

Flow at distance from wall with variable shear stress across 
the tube . - A result that avoids the assumption of uniformity of 
shear stress across the tubs was derived by von Kannan substantially 
as follows: The relation for the variation of shear stress with 

radius for fully developed flow is obtained by equating the shear 
forces to the pressure forces acting on a cylinder of fluid of 
arbitrary radius and differential length (reference 5). This rela- 
tion is 


T = — T q = (i - X) T 
r n 0 r„ 0 


or 


T = (1 - £— ) T 
+ U 

r o 


( 16 ) 
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where 



(17) 


c 

t 


On auhstitutlon of equations (7) to (10) and (16), equation (3) 
■becomes, in dimensionless form, 




( 18 ) 


where the negative sign was selected on taking the square root in 
order to make K positive. 

The first integration gives 


«= - 2K V^Vr 0 + - y + + C x 

du + 


As previously stated, the velocity gradient at the wall is very 
large compared with that at a distance from the wall and dy + /du + 
can be considered to be zero at the wall (that is, at y + = 0). 
Hence 

C x « 2Kr 0 + 

and 

&L = - 2* A Jr^ / \|r Q + - y + + 2itr 0 + (19) 

du + 


I 
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By substitution of J? = Tq* - y + , equation (IS) can be integrated 
to give 


vr 


1 

K 



zL + 


log. 


0 



( 20 ) 


Equation (20), which is substantially the equation obtained by 
von Karman, relates u + to y + for various values of r Q + . Taking 
the variation of shear stress into account introduces the additional 
dimensionless parameter r^*, which will be called the tube-radius 
parameter. The equations developed by von Karman for the flow 
distant from the wall will be extended for compressible flow by a 
method subsequently developed herein. 


Equations for Compressible Flow 

In the derivation of the compressible-flow equations, uniformity 
of shear stress across the tube. is assumed and the error introduced 
is discussed in the section "Effect of Variable Shear Stress." 

Flow in vicinity of wall . - For compressible flow in a tube, 
the static temperature and thus the density and the viscosity vary 
across the tube. For flow without external heat transfer, the total 
temperature at any point across the tube is practically the same as 
the wall temperature, so that the static temperature is given by 


t 



( 21 ) 


This relation is exact for a Prandtl number of 1. 

By use of the perfect gas law, the density o may be expressed 
as 


' 0 


_£_ JE 

gRt 7 u 2 \ 

88 (*0 - 



u" 


1 


2gJc p T 0 


f *.'*'*-‘1 % ’4iFBfr 


t" t 
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vhere the static pressure p has been assumed uniform 
tube. 


across the 


v v Viscosity is a function of temperature and can be represented 
by an equation of the form represented 


4 = Kt c 


vhere K is some constant. Then 


M = K T n - _H = KT d > l - u 2 

\ ° 2gJc p J 0 ^ 2g J Cp T 0 


4 = H 0 1 


2gJc p T 0 




4 0 1 - 


■~5 — n2 « y 


2gJc p T 0 
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oj 

c\) 

V »J 


Substituting values from equations (7) to (9) into equation (25) in 
order to convert it to a dimensionless form, gives 


s!h!z! + (i - au +2 ) d in! = i 

i-au +2 d y + ' d y + 


where 


a 


T 0 

2gJc p T Q p 0 


(27) 


The dimensionless parameter a vill be called the compressibility 
parameter. After rearrangement, equation (26) becomes 


<kl . y + = fl - 

. + .2 V 

d u 1 - au + 



(28) 


which is a first-order linear differential equation with variable 
coefficients. Solving equation (28) and setting the constant of 
integration equal to zero, as was done in the solution of equa- 
tion (11), gives 



^ 0 


For compressible flow in the vicinity of the wall, equation (29) 
gives the relation between u + and y + for various values of a. 


Flow at a distance from wall . - With the assumption of uniform 
shear stress, equation (3) can be written for compressible flow as 
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I ^- u 

.2 Uy/ J 


1 u^ / d2u \^ 

“ 2gJc p T 0 \^2/ 

When written in dimensionless form, this equation "becomes 


dV 


II - au 


+ 2 w* 


One integration of equation (30) can be made by substituting 
v = du + /dy + and dy + = du + /v and then separating the variables. 
After one integration, equation (30) becomes 


log (C, — ^ = - -pn sin" 1 (Wau + ) 
6 ^ 1 dy V A/a 


+ - 77 = sin- 1 (Wau + ) 

i _ ^ A/a 


rp= sin" 1 (A/au + ) 

+ _ o.-'Na *„+ 


dy * C-ie 


If — sin" 1 ( J Vau + ) is set equal to Z, equation (32) can be 
A/a 

integrated to yield 
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i*. 1 


K 2 +a 


-£= sin"-^ (Vau + ) 

N<x 


! J(kC 


au + ^ + au + J + C 


■~= Bin ' 1 (W) 

By substituting dy /du from equation (31) for C-^ e a 

In the immediately preceding equation and letting dy + /du + equal 

zero when y + and u + equal zero, G equals zero and thus 


+ Ci 

y *= 


K 2 +a 


e 


-£= sin -1 (Vau + ) 
AJa 


VTT 


+*• + 
au + au J 


(33) 


For compressible flow at a distance from the wall, equation (33) 
shows the relation between u + and. y + for various values of a. 


Summary of Velocity-Distribution Equations 

The Important equations that have been derived are assembled 
here for convenience: 

For incompressible flow close to the wall, 


1 

V2rt 


,+ 1 


/-mi 1 


0 


+ %2 


-(nu+) 
e 2 d(nu + ) 


-(nu+) 


+ \2 


( 12 ) 


1 e 2 


V2x 


t 

-( nu ) - 

_1_ e 2 is the normal error function of nu + . 
V2rt 
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For incompressible flow at a distance from the wall, with shear 
stress uniform across the tube, 

u+ c 7 log e y + + C (15) 

For Incompressible flow at a distance from the wall with shear 
stress variable across the tube, 


u + = i 
K 


£— + iog A li -Ji - zL 


+ c 


For compressible flow close to the wall. 


y + = 


1 - au 


-r? 

j\ 2a 


pu 1 


1 - au 


+ 2 V a 


i 2 ,, 


du 


( 20 ) 


(29) 


For compressible flow at a distance from the wall with shear stress 
uniform across the tube, 


K^+a 


M 


sin - ^ (/Vau + ) 



(33) 


APPARATUS 


A schematic diagram of the experimental apparatus used is showr. 
in figure 1. Air at a pressure of about 40 pounds per square inch 
gage and at approximately room temperature flowed through two con- 
trol valves in parallel, then through a filter, an orifice, a calming 
tank, a test section, and into the atmosphere. 
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Test section . - The te3t section consisted, of a smooth Inconel 
tube having an inside diameter of 0.87 inch, an outside diameter of 
1.0 inch, and a length of 87 inches. The calming tank was so con- 
structed that either a rounded or a right-angle-edge entrance to the 
tube could be used; the rounded entrance consisted of a standard 
A.S.M.E. long-radius nozzle. Static pressures were measured through 
0.03- inch holes drilled in the tube at the positions shown in 
figure 1. Air total temperatures were measured by thermocouples 
located at the inlet and outlet of the tube; these measurements 
provided a check on whether heat transfer took place through the 
tube wall. 

Total-pressure measurements . - Openings fcr taking total- 
pressure measurements across the tube were located as shown in 
figure 1. Holes having 0.15-inch diameters through which a total- 
pressure probe entered the tube were drilled in the tube wall. 

These holes were located at right angles to the static-pressure 
taps. A probe actuator to move the probe across the tube and to 
measure its distance into the tube was fitted into a short length 
of tubing at each opening. The location of the total-pressure 
probe with respect tc the opening in the tube is shown in the insert 
in figure 1. The total-pressure probe used for the measurements had 
a 0.016- inch -diameter opening and a 0.002-inch wall at the tip. The 
probe was flattened out for some of the runs so that the width of 
the opening was 0.005 inch. The probe was made so that the tip 
just cleared the edge of the 0.15- inch hole in the test section. 

The total projected area of the probe in the direction of flow 
with the tip at the center of the tube was about 1-i percent of 
the area of the tube, but the effective blocking area at the tip 
was considerably less because the main portion of the probe was 
downstream of the tip. 


METHODS 

Procedure 

In order to establish the applicability of the equations derived 
for fully developed flow, velocity distributions at various dis- 
tances from the tube entrance were first determined. Measurements 
were made at Reynolds numbers of approximately 40,000, 160,000, and 
580,000 with both rounded and right-angle-edge entrances and with 
corresponding tube-exit Mach numbers up to 1. The flow rate in 
each case was obtained by adjusting the inlet pressure. 
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A more extensive series of investigations vas then conducted to 
determine velocity distributions for fully developed flow including 
distributions for the region close to the wall. Runs were made at 
about 20 different Reynolds numbers between 10,^00 and 200,000 and 
with Mach numbers up to 0.5 with the tctal-pi'essure probe near the 
exit of the tube where the flow was practically fully developed. 

Runs were made with both the rounded and the right-angle -edge 
entrances. Total- pres sure readings were taken at points between 
the wall in which the probe opening was located and the center of 
the tube. Readings near the opposite wall were inaccurate because 
of disturbance due to the probe. 

Measurements were made with both the round and the flattened 
probe tips in order to determine whether presence of the probe had 
any effect on the measured velocity distributions in the vicinity 
of the wall. Because no appreciable difference could be noticed 
between the measurements made with the two types of tip, it was 
concluded that the presence of the probe did not affect the meas- 
urements and that the aerodynamic and geometric centers of the hole 
coincide . 

Preliminary runs at high and at low flow rates shewed that the 
total temperature of the air was uniform along the length of the 
tube, indicating that no heat transfer occurred. The air total 
temperature was therefore measured only at the cutlet. 

The following quantities were measured for each run: air flow, 

static pressures at the wall, air total temperature, total pressure 
at various positions across the tube, and distance of the probe tip 
from the wall. 


Reduction of Experimental Data 

Velocities . - For low air-flow rates, incompressible-flow 
theory was used; the velocities were calculated from the equation 


P 



pu 


where 


o was found from the equation of state for perfect gases 
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and t was taken equal tc the total temperature. In this and in all 
succeeding calculations, the static pressure was assumed to be uni- 
form across the tube. 


For Mach numbers greater 
from the relation 


than 0.2, 


velocities were calculated 


P 

P 


+ 


Hi 

7 



111 

where t = T ^ 

Shear stress. - The shear stress at the wall for fully dev- 
eloped flow is related to the friction- pressure gradient by the 
equation 



The friction-pressure gradients were obtained by subtracting 
calculated momentum-pressure gradients from the measured static- 
pressure gradients along the tube; the moment urn- pres sure gradients 
were calculated from 


V mom p^A^g^ ^ 


where 


Pb 


was found from the equation of state 


P = 



> *&*> V.JUaJrSfc'fc. • ►*• w. . 
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and 


■1 +>1 + 2 


^ = 


w 2 R 2 T 

gJCpA 2 p 2 


v 2 R 2 


gJc p A 2 p 2 


The last equation was obtained from the equations of energy, con- 
tinuity, and state. The pressure and density gradients were 
graphically determined by plotting pressure and density against 
distance along the tube and drawing a tangent to the curve at the 
point in question. 

Distance from wall . - The zero reading of the probe actuator 
was found by plotting velocity against distance reading on the 
probe actuator for a number of runs and extrapolating the curves 
to zero velocity where all the curves Intersected. This extrapola- 
tion gave the probe-actuator reading with the probe tip at the wall. 
The distance of the tip from the wall for each reading could then 
be easily calculated. 


Bulk velocity. - In order to obtain ratios of local to bulk 
velocity ,u / u^ at various positions across the tube, the bulk 
velocity u^ was obtained by plotting u against cross-sectional 
area nr 2 measuring the area under the curve, and dividing by 


the total cross-sectional area of the tube 
is equivalent to solving the equation 


itr. 


This procedure 






u d(nr 2 ) 


nr 


0 


and gave more accurate values of u/u-^ than would have been 

obtained by use of orifice measurements of weight flow for the 
determination of u^, Inasmuch as errors in u due to errors 

in static-pressure measurements were also contained in u-^ and 
any systematic errors tended to cancel. 
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RESULTS AND DISCUSSION 

Variation of Velocity Distribution along Tube 

The results of the tests in which variation of the velocity 
distribution along the tube was determined are summarized in 
figures 2 and 3. Typical radial velocity distributions at various 
tube diameters from the rounded entrance of the tube are shown in 
dimensionless form in figure 2. 

In figure 3, which Is somewhat mere descriptive than figure 2, 
the variation of velocity along the tube at the center and at 
t/yq =0.9 is shown for both a rounded and a right-angle -edge 
entrance. The curves show that the development of the velocity 
distribution was more rapid for the right-angle-edge entrance than 
for the rounded entrance; with the right-angle-edge entrance fully 
developed flow was obtained after about 45 tube diameters from the 
entrance, but with the rounded entrance the distribution was still 
developing slightly at 100 tube diameters from the entrance. This 
difference In rate of development wa3 apparently caused by the vena 
contracta formed at the entrance of the tube with the right-angle - 
edge entrance, and indicated on figure 3(b) by the points close to 
the entrance. The presence of the vena contracta accelerated the 
flow at the center of the tube and thus hastened the development of 
the distribution. 

The curves also show that Reynolds number affected the dis- 
tribution near the center of the tube because u/u-^ decreases as 
Reynolds number increases. This variation Is in agreement with 
previous findings (reference 6). Difference in Reynolds number 
has, however, little or no effect on the distribution close to 
the wall. 

A significant observation can be made from figure 3 concern- 
ing the difference between the rates of development of velocity 
distribution at the center of the tube and near the wall. Figure 3 
indicates that with both entrances the final distribution is 
attained much sooner in the vicinity of the wall than at the 
center of the tube. This fact might explain why, in the present 
investigation, the static-pressure gradients along the tube caused 
by friction were only slightly affected by entrance effects. The 
static-pressure gradient is determined by the velocity gradient 
at the wall and is unaffected by the distribution in the remainder 
of the tube. 
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As will be shown in the discussion of fully developed flow in 
the section "Effect of Variable Shear Stress," the effect of com- 
pressibility on velocity distribution is slignt; although some of 
the determinations were made at high subsonic Mach numbers, it is 
therefore believed that the results are also applicable to incom- 
pressible flow. 


Velocity Distributions for Fully Developed Flow 

Correlation of experimental data . - The variation of u + with 
y + for data obtained near the exit of the test section (that is, 
for x/D = 100) where the flow is fully developed is shown on 
rectangular coordinates in figure 4(a). The data obtained near 
the wall are plotted to two y + Beales. The data are plotted 
semilogarithmically in figure 4(b). Data for flew close to the 
vail are shewn for only the lew flow rates because at high flow 
rates the severe velocity gradients and the presence of the hole 
In the tube wall make the accuracy of the measurements doubtful. 

The data shown were taken with both rounded and right-angle-edge 
entrances, but the velocity distributions with the two types of 
entrance were the same within the error of the measurements. 

Comparison of the data in figure 4(b) with those of Nikuradse 
and of Reichardt (reference 3) shows close agreement for all values 
of y + up to about 600. For higher values of y + , the corres- 
ponding values of u + are slightly higher than those obtained 
by Nikuradse (reference 1, p. 242); the maximum deviation, however, 
is only about 5 percent. 

A method of using figure 4 to obtain the velocity distribution 
for a particular tube when the flow rate is given is indicated in 
the discussion of figure 8. 

Incompressible-flow equations . - The curve corresponding to 
equation (12) for incompressible flow near the tube wall is in- 
cluded in figure 4 and is in good agreement with the experimental 
results for values of y + from 0 to 23. The value of the con- 
stant n in the equation is 0.109, as determined from the 
experimental data. 

An important property of equation (12) is that for small 
values of y + , u + and y + are approximately equal, that is, 
the flow predicted by the equation is nearly laminar. This approach 
to equality of u + and y + is due to the fact that as the wall 
•a approached the eddy diffusivity € becomes very small and is 


K 
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zero at the vail, (equation (5)). The accuracy with which equa- 
tion (12) predicts the experimental data for y + from 0 to 26 
obviates the necessity of assuming the existence of a separate 
layer near the wall that la piirely laminar, but does not eliminate 
the possibility of its existence, as has usually been done in pre- 
vious investigations. A single equation has therefore been obtained 
that, for adiabatic incompressible flow, adequately represents the 
two regions, which are commonly called the laminar layer and the 
buffer layer. The buffer layer has previously been represented only 
by empirical equations. 

The agreement of equation (12) with the data for values of 
y + from 0 to 23 does not eliminate the possibility of the exist- 
ence of a very thin layer that is purely laminar, for example, for 
the region 0 < y + < 3. It is possible that a finite layer exists 
in which adjacency of the wall completely damps out turbulence. 

The thickness of the layer then corresponds to some critical wall 
Reynolds number ou e y Q /p where u 8 is the velocity at the edge 
of the layer and y 9 is the thickness of the layer. Velocity 
distributions for 0 < y + < 3 given by equation (12) are prac- 
tically laminar, however, so that it makes little difference, for 
calculating velocity distributions, whether the layer is taken 
Into consideration. The only case in which presence of the layer 
may become important Is that of heat transfer in fluids having 
high Prandtl number, where the turbulence predicted by equation (5), 
though it may be slight, is important because of the small amount 
of heat transferred by conduction. 

The agreement of equation (12) with the data, together with 
the discussion preceding the derivation of the equation, indicate 
that in the region close to the wall the mechanism of turbulent 
transfer of momentum can be considered affected mainly by quanti- 
ties that are determined relative to the wall; that is, by the 
distance of the point from the wall and by the velocity at the 
point relative to the wall. As was shown in the discussion pre- 
ceding the derivation of equation (12), the velocity distribution 
about the point is known to a first approximation when the two 
quantities u and y are known at the point. 

The general form of the equation that is usually employed to 
represent the turbulent regime was obtained by von Karman and is, 
as shown in the analysis section, 


u + = I lcg e y + + C 

ft 


( 15 ) 
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The values of the constants C and K, obtained from 
Nikuradse's data (roference 1), ore 5.5 and C.40, respectively. 
Corresponding values that represent the present data somewhat 
better are C = 3.S and K = 0.36. The curve representing 
equation (15) using these values is plotted in figure 4, in which 
good agreement with the experimental data is indicated for y + = 0 
to y + - 26. The curves corresponding to the equations for flow 
near the wall (equation (12)) and flew at a distance , from the wall 
(equation (15)) cross at y + = 26. 

The agreement of equation (15) with the data, together with 
the discussion preceding the derivation of the equation, indicate 
that for a region distant from the wall the mechanism of turbulent 
transfer of momentum can be considered dependent mainly on the flow 
conditions in the vicinity of the point considered: that is, on the 

velocities in the vicinity of the point relative to the velocity at 

the point and not on the position of the point in the tube (the 

distance from the wall) or on the velocity relative to the wall. 

The curves representing equations (12) and (15) have slopes 
that are not quite equal at their Intersection at y + = 25: this 
disparity would, however, be expected because the two equations 
were derived with the assumption that the turbulence mechanism in 
the two regions is dominated by different factors: hence an abrupt 
change in turbulence mechanism at the intersection is implied in 
the equations. Actually, there is probably a gradual change that 
could not be investigated by the simplified methods used herein. 
Inasmuch as the actual error in the vicinity of the intersection 
is insignificant, the present treatment is considered adequate for 
adiabatic flow. 

The value for K (in equation (15)), which is known as the 
Karman constant, was checked by calculating friction factors and 
Reynolds numbers and plotting l/Y4f against ReV4f as 3hown 
in figure 5. The line drawn through the data corresponds to the 
Karman relation between friction factor and Reynolds number, which 
is 


_A_ = Co + log (Se Vif ) 

Vif k •fa 


This equation is derived In reference 2 directly from the 
equation for velocity distribution. The value of K vas 0.36 
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as before , and C 2 Is found from the data to be - 1.84. Both 
velocity-distribution and friction- factor data therefore indicate 
that a value for K of 0.36 is reasonable, at least for the 
present tube. 

The variation of the ratio of velocity at the center of the tube 
to bulk velocity with Reynolds number, as indicated in figures 3(a) 
and (b), can be explained by the plots of the equations in fig- 
ure 4(a). As y+ increases the curve becomes flatter. Increasing 
values of Reynolds number correspond to increasing values of y + 
in the central portion of the tube, so that for high values of 
Reynolds number the velocity profile in the central portion of the 
tube becomes flat; thus the ratio cf velocity at the center of the 
tube to bulk velocity becomes less than for low Reynolds numbers. 

Effect of variable shear stress . - Neither u + nor y + are 
functions of the tube radius. Neglecting the radius in the correla- 
tion is equivalent to assuming uniform shear stress across the tube. 
(See equations (15) and (20) in the section "Analysis . " ) The 
results can be correlated by assuming uniform shear stress across 
the tube because the greatest rate of velocity change with dis- 
tance from the wall occurs near the wall where the change in shear 
stress is very small, A comparison between equation (20), which 
takes into account the variation in shear stress across the tube, 
and equation (15), which assumes uniform shear stress, is shown in 
figure 6. The constant C is so determined for each value of 
Tq+ that u+ = 13 when y + = 28. These values were selected 
in order to make the mean deviation of equation (2C) from equa- 
tion (15) a minimum. The maximum difference between the values 
of u + determined by the two equations is about 2^ percent, which 

is the same as the scatter of the experimental data points. Equa- 
tion (15) therefore gives an accuracy comparable to that obtained 
in most flow measurements so that consideration of the variation 
in shear stress across the tube does not seem necessary. 

Equations for compressible flow . - The equations for com- 
pressible flow are presented in figure 7. For graphing equation (29), 
the value of the integral was found by plotting the integrand against 
u+ and planimetering the area under the curve. The constant n 
was again set equal to 0.109 and d for air was found from vis- 
cosity data to have an average value of 0.684 for temperatures 
between 0° and 200 0° F. For plotting equation (33), was deter- 

mined for each value of a from the value of u + at y + = 26 
found from equation (29), From the definitions of M, T, u + . 
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and a, It can be shown that the relation between u + , M, 7 , 
and a Is 


“(i^r-i) * *) 

For M = 1, 7 = 1.4, and a = 0.00025, u + Is found to he 25.8. 

This point is marked on the curve. 

As a Increases, u + decreases for constant values of y + 
(fig. 7). For local Mach numbers up to 1, hoyever, the deviation 
is slight and, in general, is not more than 2^ percent from the 
value of u + given by the incompressible-flow equation. A single 
line is used to represent equation (29) because the compressibility 
effect Is so small that it cannot be seen when the equation is 
plotted to the scale used in this figure. 

It therefore appears that with respect to compressibility 
effects and tube radius or shear-stress variation, the simplified 
incompressible-flow equations (12) and (15) give an accuracy com- 
parable to that of flow measurements. The fluid properties used 
in u + and y + are evaluated at the wall or total temperature. 


Flow Rates 

The flow rates corresponding to various pressure gradients 
along the tube can be obtained from the velocity-distribution 
equations by a graphical integration. For this Integration, 
(p/pq) u + was plotted against (r 0 + - y+) 2 , where 
p/Pq = l/(l-au + j. It can easily be shown that the area under 
this curve is w (A/t q/ 0 q pQ^jtgUQ 2 . If this dimensionless 
group is divided by Tq + and multiplied by it in order to 
eliminate T^, there results w/(gpQrQ). This parameter is 
plotted against r 0 + for various values of a in figure 8. 
Figure 8 gives the flow rate to be expected for a given shear 
stress or friction- pressure gradient. The data of figure 8 also 
provide a means for obtaining the velocity distribution in a 
particular tube from the generalized velocity distributions in 
figures 4 and 7 when the flow rate and the fluid properties are 
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known. Measured, flow rates for conditions where compressibility 
effects are small are also plotted, in figure 8 and. agree closely 
with the curve for a=0. The curves indicate that at a Mach number 
of 1, compressibility effects increase the flow-rate parameter 
w/(gg n r n ) by about 12 percent. Several data points corresponding 
to higher values of a are also plotted in the figure and the trends 
corresponding to increasing values of a appear to be similar to 
the predicted trends, although no definite conclusion can be drawn 
from the limited range of values of a shown. Data corresponding 
to higher values of a are not plotted because reliable measure- • 
ments of the severe pressure gradients involved could not be 
obtained. 


CONCLUDING REMARKS 

The results obtained in this Investigation should be applicable 
to any gas to which the perfect gas law applies and for which the 
Prandtl number is close to 1. The value of the exponent d for 
viscosity variation with temperature was obtained specifically for 
air, but it occurs only In the equation for compressible flow close 
to the wall where compressibility effects are negligible. The 
equations and the curves for Incompressible flew should, of course, 
also be applicable to liquids, inasmuch as the fluid properties 
that determine the flow phenomena are common to liquids and gases. 


SUMMARY OF RESULTS 

The following results were obtained from the analytical and 
experimental investigation of the adiabatic flow of air through a 
smooth tube having an inside diameter of 0.87 inch and a length of 
87 inches : 

1. The length of tube required for obtaining fully developed 
flow was greater with a rounded entrance than with a right-ang.le- 
edge entrance. With a rounded entrance, the flow at the axis of 
the tube was still developing slightly at 100 diameters from the 
entrance. For both entrances, however, the flow close to the wall 
developed In a much shorter distance than did the flow in the center 
of the tube. The flow close to the wall determines the shear stress 
or pressure gradient along the tube, so that the effect of entrance 
on these factors is Blight except for very short tubes. 
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2. A good correlation of velocity-distribution data for fully 
developed flow was obtained by using the well-known dimensionless 
velocity and distance parameters. The data agreed closely with 
these of Nikuradse and other investigators. 

3. The equation derived for adiabatic incompressible flow close 
to a wall represented well the relation between the velocity and 
distance parameters found experimentally for the two regions 

that have generally been called the laminar layer and the buffer 
layer , 

4. The analysis and experimental investigations indicated that 
the effect cf variable shear stress on velocity distributions is 
slight; the maximum variation in the velocity parameter caused by 
this effect was approximately 2^ percent. 

5. The compressible-flow equations showed that the effects 
of compressibility on velocity distributions are small for Mach 
numbers up to 1: the maximum variation in the velocity parameter 
caused by compressibility effects va3 approximately 2^ percent. 

6. The simplified incompreeeible-ficv equations derived on 
the assumption of uniformity of shear stress across the tube pre- 
dicted velocity distributions in smooth tubes for Mach numbers up 
to 1 with an accuracy comparable to that of flow and pressure 
measurements . 

7. Graphical integration of the velocity-distribution equa- 
tions gave flew rates that agreed closely with flew rates from 
orifice measurements. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 9, 1950. 
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APPENDIX - SYMBOLS 


The follcving symbols ere U3ed in the report: 

internal cross-sectional area of tube, sq ft 


C > C 1j C 2 


constants of integration 

specific heat of fluid at constant pressure, Btu/(lb)(°F) 
inside diameter of tube, ft 

exponent that describes variation of viscosity of fluid 
vith temperature 

acceleration due to gravity, 32.2 ft/sec^ 
mechanical equivalent of heat, 778 ft-lb/Btu 
constant 

total pressure, lb/sq ft absolute 
static pressure, lb/sq ft absolute 
perfect gas constant, ft-lb/(lb) (°R) 
radius, distance from tube center, ft 
inside tube radiuB, ft 
total temperature, °R 
absolute wall total temperature, °R 
absolute static temperature, °R 

bulk or average static temperature of fluid at cross 
section of tube, °P. 

velocity parallel to axis of tube, ft/sec 

bulk or average velocity at cross section of tube, ft/sec 

velocity at center of tube, ft/sec 
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fluid-flov rate, lb/sec 

axial distance from tube entrance, ft 

distance from tube vail, ft 

ratio of specific heat at constant pressure to specific 
heat at constant volume 

coefficient of eddy diffusivity, sq ft/sec 

Karman constant 

absolute viscosity of fluid, lb-sec/eq ft 
absolute viscosity of fluid at vail, lb-sec/sq ft 
mas 8 density, lb-sec 2 /ft 4 

bulk or average density at cross section of tube, 
lb-sec 2 /ft 4 

mass density of fluid at vail, lb-sec 2 /ft 4 
shear stress in fluid, lb/sq ft 
shear stress produced by turbulence, lb/sq ft 
shear stress produced by viscosity, lb/sq ft 
shear stress in fluid at vail, lb/sq ft 


Subscripts : 


on friction pressure gradient 
on momentum pressure gradient 


Dimensionless parameters: 


compressibility parameter, 


2gJc p T 0 P 0 
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friction factor, - 


E (&)fr 

2 (>b u b 2 


M 


Mack number, — . 

Vt sRt 


Re 


Reynolds number. 


Pb u b D 

^0 



tube-radius parameter. 


vy pp 

V *>0 


r 


0 


u + velocity -parameter. — H 

vy^ 


vall-distance parameter. 


V T o/ °o 

y y 


y 
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Figure 1. - Apparatus for measuring velocity distributions 
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Inside tube radius' rg 

Figure 2. - Typical velocity distributions for flow through tube 

with rounded entrance at various distances from entrance. Reynolds 
number, 46,000. 
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Figure 3* - Development of velocity distributions for flow through tube 













Figure 4. - Concluded. Generalized velocity distribution for fully developed flow in smooth tubes 




Figure 5. - Determination of the Karman constant X from friction-factor data 






Comparison of equation derived by assumption of uniform shear stress with equation accounting 
for variation of shear stress across tube* 
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Figure 7* - Velocity-distribution equations for compressible flow in smooth tubes 





